Low-temperature properties of classical, geometrically frustrated antiferromagnets 
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We study the ground-state and low-energy properties of classical vector spin models with nearest- 
neighbour antiferromagnetic interactions on a class of geometrically frustrated lattices which includes 
the kagome and pyrochlore lattices. We explore the behaviour of these magnets that results from 
their large ground-state degeneracies, emphasising universal features and systematic differences be- 
tween individual models. We investigate the circumstances under which thermal fluctuations select 
a particular subset of the ground states, and find that this happens only for the models with the 
smallest ground-state degeneracies. For the pyrochlore magnets, we give an explicit construction 
of all ground states, and show that they are not separated by internal energy barriers. We study 
the precessional spin dynamics of the Heisenberg pyrochlore antiferromagnet. There is no freezing 
transition or selection of preferred states. Instead, the relaxation time at low temperature, T, is of 
order ft/fesT. We argue that this behaviour can also be expected in some other systems, including 
the Heisenberg model for the compound SrCrsGa^OiQ. 



PACS numbers: 75.10.Hk, 75.40. Mg, 75.40.Gb 
I. INTRODUCTION 

Experimental and theoretical studies in recent years 
have found that geometrically frustrated antiferromag- 
nets display properties quite unlike those of other mag- 
netic systems.El These materials have magnetic ions lo- 
cated on lattices of site-sharing frustrated units - usu- 
ally triangles or tetrahedra. One of the best-studied sys- 
tems in this class is the layered compound SrCrgGa^Oig 
(SCGO) .EtEI Attention has focussed on the fact that the 
majority of its magnetic CY 3+ -ions reside on the sites 
kagome lattices, although the full structure is more com- 
plex. Following the interest in kagome magnets gen- 
erated by studies of SCGO, a great deal of attention 
has been devoted to the oxide and fluoride pyrochlore 
magnets, in which the magnetic ions form a lattice of 
corner-sharing -tetpahedra as depicted in Fig. LTL r-Neu- 
tron scatteringoBii3 and muon spin relaxation&EJ ex- 
periments on SCGO and the pyrochlores have detected 
only short-range magnetic correlations and a slowing- 
down of fluctuations at low temperatures .□ More gen- 
erally, it is apparently a characteristic property of ge- 
ometrically frustrated magnets that they do not order 
at the temperature expected from the magnitude of the 
Curie- Weiss constant, |0cw|- Instead they remain in the 
paramagnetic phase to a much lower tempeca,ture with - 
typically - spin freezing at Tp < |0cw| •circa 

A detailed understanding of the origin of such generic 
features has been slow to emerge. Moreover, there has 
been little work to explain systematic differences between 
individual examples of these magnetic systems. For in- 
stance, whereas in the Heisenberg kagome antiferromag- 
net thermal flipteations give rise to entropic growid-. 
state selection,c-iEa known as order by disorder ,E3~E2l 
this phenomenon appears to be absent for some related 
systemsH'El The reason for the difference is unclear, as 
are the general conditions under which such selection 
should be expected for geometrically frustrated magnets. 
In this context, it is interesting to ask whether SCGO 



inherits its properties from those of the kagome Heisen- 
berg antiferromagnet, or whether its behaviour is closer 
to that of the pyrochlore antiferromagnet, since an al- 
ternative and more complete description of the structure 
is to regard a layer of SCGO as a slab cut from the 
pyrochlore lattice, consisting of three consecutive [111] 
lattice planes. 




FIG. 1. The pyrochlore lattice 

A further area for investigation, in addition to the sta- 
tistical mechanics of geometrically frustrated antiferro- 
magnets, is their low-temperature cbmamics, which has 
so far received only limited attention.Eja Dynamical cor- 
relations are likely to be profoundly influenced by the 
large ground state degeneracy of these systems, and con- 
stitute one of their most interesting aspects. 

In an attempt to extend understanding of these prob- 
lems, we have studied the low-temperature properties of 
the classical Heisenberg model with nearest-neighbour in- 
teractions on a class of geometrically frustrated lattices. 
This description neglects various-additional features - 
such as anisotropyj2jti-disprder,Ej dipolaio or jsAT|thft&. 
neighbour interactions ,e3c3 and quantum effect s£3'E3 E2l 
- which can play an important role in real materials, par- 
ticularly near and below Tp. However, it may provide a 
good treatment for the temperature window T < Tp <C 
|@CvH) ano - simplicity should make it well-suited for 
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capturing the generic features of these systems, as well as 
providing a basis for future investigations incorporating 
additional interactions or quantum fluctuations. 

Our work, parts of which have been described in 
Refs. |5l| and [32| concentrates on the pyrochlore antifer- 
romagnet, but we address several questions in a more 
general context. We start by analysing the origin and 
extent of the ground-state degeneracy of geometrically 
frustrated magnets (Sect. ||). We discuss the nature of 
the ground-state manifold of pyrochlore antiferromag- 
nets with rt-component spins (Sect. III). We give an 
explicit construction of all ground states of these mag- 
nets, and show that they are not separated by energy 
barriers. Although typical ground states are disordered, 
we show that certain correlations remain, which give 
rise to distinctive features in magnetic neutron scatter- 
ing. We study, both analytically (Sect. IV) and numer- 
ically (Sect, [v]), the existence of order by disorder for 
a general class of geometrically frustrated antiferromag- 
nets and find that it occurs only for magnets with small 
ground-state degeneracies. In particular, it is absent 
from the Heisenberg pyrochlore magnet, which therefore 
has neither internal energy nor large free energy barri- 
ers separating different ground states. Because of this, 
the system is not trapped near a particular state at low 
temperatures. Our study of the precessional dynamics 
at low temperatures and for long times (Sect. VI) re- 
veals that the decay of the autocorrelation function is 
exponential in time, t, with a timescale inversely pro- 
portional to the temperature and independent of the ex- 
change energy: (Sj(O).Sj(i)) = exp(— cksTt/h), where c 
is 0(1). In agreement with Reimers' earlier Monte Carlo 
simulations ,E3 we find that the spin- freezing transition 
observed experimentally does not happen in the simple 
Heisenberg model we consider. ,-We discuss recent exper- 
iments on pyrochlore mag netsMl and SCGOVM in the 
light of these results. 

Since spin correlations are short-ranged in both space 
and time, the Heisenberg pyrochlore antiferromagnet can 
be labelled a classical spin liquid or, following Villain,E3 
a cooperative paramagnet. 



II. THE HEISENBERG SPIN HAMILTONIAN ON 
GEOMETRICALLY FRUSTRATED LATTICES 

Consider n-component classical spins, Si, with |Si| = 
1, arranged in corner-sharing units of q sites. Each spin 
is coupled antiferromagnetically with its q — 1 neighbours 
in each unit, so that the Hamiltonian is 



H = J^ S ? -S^^|L q | 2 - J -Nq. 



(2.1) 



<i,3> 



Here, J is the exchange constant and L a is total spin in 
unit a. The sum on runs over all neighbouring pairs 
and the sum on a runs over the N units making up the 
system. 



Note that our motivation for considering n-component 
spins is to shed light on the systematics of geometrically 
frustrated antiferromagnets. Because of this, we take the 
rt-component spin-space to be the same at each site. Of 
course, the case n = 2 can also arise physically in a 
Heisenberg system with easy-plane anisotropy: in this 
event, which has been studied in Ref. [54|, the easy planes 
are orientated differently at different sites, in accordance 
with the local symmetry axes. 

An instructive way of thinking about the strength of 
the geometric frustration is to consider the extra ground- 
state degeneracy which it gives rise to, in addition to 
that stemming from the symmetry of the Hamiltonian. 
It is this extra degeneracy which lies behind many of the 
physical properties peculiar to geometrically frustrated 
systems. To determine the number, D, of degrees of free- 
dom in th e gr ound state, we use a Maxwellian counting 
argument ,e3Ej and evaluate Dm = F — K, the difference 
between the total number, F, of degrees of freedom in 
the system, and the number, K, of constraints that must 
be imposed to restrict the system to its ground states. In 
general, as discussed below, Dm ^ D, b ut for pyrochlore 
antiferromagnets we argue in Sect. [II B that Dm/D — > 1 
as N — > oo. 

To evaluate K, note that, from Eq. 2.1, a configuration 
is a ground state provided L Q = for each unit sepa- 
rately. This imposes Nn constraints. To find F, we start 
from the fact that the number of degrees of freedom is 
simply n—1 per spin. Expressed in terms of the number, 
N, of units, F depends on their geometric arrangement. 
For corner-sharing units of q spins, F — Nq(n — l)/2. 
Alternative arrangements generally result in smaller val- 
ues of F/N and in ground states that are not extensively 
degenerate. For example, if bonds are shared between 
units - as in the triangular and face-centred cubic lat- 
tices for q — 3 and q — 4 respectively - F is lower than 
if only sites are shared - as in the kagome and the py- 
rochlore lattices - since each spin belongs only to b = 2 
units in the latter case but to more (b = 6 and 8, re- 
spectively) in the former. In the general case, we obtain 
F = Nq{n- l)/b. Hence, D M /N = [q(n - l)/b - n). 
Dm grows with q and, for q > b, with n. In order to 
obtain Dm > 0, we require q > b, which is the case only 
for corner-sharing arrangements. The physically realis- 
able example for which Dm is maximal is that for which 
q and n are both maximal: Heisenberg spins (n = 3) 
on the pyrochlore lattice (q = 4) represent the only sim- 
ple system for which Dm is positive and extensive. It 
is partly for this reason that the pyrochlore Heisenberg 
antiferromagnet is particularly interesting. 

This counting argument can go wrong in two ways. 
Firstly, the K constraints may not be independent, as 
happens for Heisenberg spins on the kagome lattice, 
where Dm = but an extensive ground-state degener- 
acy nonetheless arises. Secondly, for some lattices there 
may be no spin configurations that satisfy the conditions 
L Q = for all a. 

Many of the results presented in this paper do not de- 
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pend on the details of the lattice under consideration but 
rather on the size q of the corner-sharing units. We find 
it useful to consider, in addition to the pyrochlore laU. 
tice, the two-dimensional square lattice with crossingst3 
(Fig. |^), which is not known to occur in nature but is 
easy to visualise. Like the pyrochlore lattice, from which 
it can be obtained by a projection in a (001) direction, it 
has q = 4 and, with Heisenberg spins, Dm — N. 




FIG. 2. The square lattice with crossings. Both solid and 
broken lines denote exchange interactions. Spins reside on the 
intersections of the solid lines. 



posing extra constraints on the spin arrangement, in ad- 
dition to the requirement that it be a ground sate. A sim- 
ple example is the set of four-sublattice states, in which 
the four spins of each unit cell are arranged to be oriented 
the same way everywhere. Any four-spin arrangement 
that is a g round state for the single tetrahedron (see sec- 
tion III A ) yields a ground state for the entire system by 



periodic repetition. Of these states, those with two spins 
parallel and two antiparallel to_a given axis are the sim- 
plest conceivable ones. VillairH has described a larger 
ground-state submanifold for the Heisenberg model, in 
which the spins of each tetrahedron form two antipar- 
allel pairs. It turns out that for XY mod el all ground 
states are of this kind, as described in Sect. HIE. 



In the following, we present complete constructions of 
the ground states for classical antiferromagnets with n- 
component spins on the pyrochlore lattice. We also show 
that the ground-state manifold is connected. We then 
examine the conequences of spin correlations in typical 
ground states for elastic neutron scattering. We conclude 
this section with a discussion of the nature of the ground- 
state degrees of freedom in such magnets. 



Also, more complicated corner-sharing arrangements 
of frustrated units are possible. Of particular experimen- 
tal importance, as mentioned above, is the combination 
of triangles and tetrahedra found in SCGO, which is de- 
picted in Fig. S. 




FIG. 3. Projection of the sites of the magnetic Cr 3 + ions m 
one layer of SCGO. The ions occupy sites of the kagome lat- 
tices in the top (solid lines) and bottom (dashed lines) planes. 
In the middle plane, the ions are located on a triangular lattice 
(indicated by circles) . All lines denote exchange interactions; 
there are also interactions between a spin in the middle plane 
and those in the two triangles which enclose it. 



III. THE GROUND STATES OF CLASSICAL 
ANTIFERROMAGNETS ON THE PYROCHLORE 
LATTICE 

It has been realised for a long time that antiferro- 
magnets on the py . rq chlore lattice have a vast ground- 
state degeneracyj£lll^3 but no explicit construction of the 
ground states has as yet been available. The nature 
of various submanifolds of the ground-state manifold is 
however known. These submanifolds are defined by im- 



A. The single tetrahedron 

The ground states of a single tetrahedron are those 
states in which the sum, L, of the four spin vectors has 
the value L = 0. In such a configuration, any two spins 
enclose the same angle as the other two. For Heisenberg 
spins, these configurations can be parametrised by two 
coordinates (e.g. a and <p m Fig. ||). The crucial feature 
is that, for any fixed a except a — 0, one can choose 
4> independently. In the special case, a = 0, if spins 1 
and 4 are antiparallel, there are two degrees of freedom 
associated with the remaining two spins, while if spins 1 
and 4 are parallel, there is no remaining freedom. These 
exceptional states, in which all spins of a tetrahedron 
are collinear, can play a central role in determining the 
thermodynamics of the system because they are favou red 



by thermal fluctuations, as discussed in section IV A 




a 



FIG. 4. Four spins of equal length with vector sum 
L = 0. Spins 1 and 4 lie in the plane of the paper, spins 
2 and 3 need not. 

For the XY antiferromagnet, there is only one con- 
tinuous degree of freedom, a, since if spins are coplanar, 
4> = 0, 7r. Ground states are therefore the configurations 
with two pairs of antiparallel spins. 
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B. The Heisenberg antiferromagnet 



2. The connectedness of the ground-state manifold 



1. The construction of the ground states 

We give in this subsection a stepwise procedure for con- 
structing any ground state of the Heisenberg pyrochlorc 
antiferromagnet, from which the number of ground-state 
degrees of freedom, D, can be determined directly. We 
also consider a similar proceedure for the square lattice 
with crossings since it is essentially the same but easier 
to explain and visualise. The idea in both cases is that 
the ground state can be built up by choosing the orienta- 
tions of spins on successive layers (planes or lines) of the 
lattice, in a way that requires no adjustments of spins in 
planes or lines already visited. We consider systems with 
open boundary conditions: in the context of this section, 
periodic boundary conditions appear to introduce signif- 
icant additional mathematical difficulties. 

We define a layer, for the square lattice with crossings, 
to be a [10] plane (Fig. ||), and for the pyrochlore lattice 
to be a [100] plane. In both cases, a layer contains the 
spins lying on equivalent edges of squares or tetrahedra - 
referred to as units from hereon - which are next-nearest, 
but not nearest, neighbours, to other units with spins in 
the layer. The spins of each unit are shared between two 
adjacent layers. Conversely, each spin belongs to a unit 
extending above and one extending below the layer. 

First, choose the orientations of the spins on the lowest 
layer of the lattice. This amounts to choosing a value for 
a in each unit with spins on the bottom layer. There 
are no restrictions on how to do this. Next, consider the 
adjacent layer: when choosing the orientation of spins on 
that layer, one has to satisfy the ground-state condition. 
For each unit, this leaves one degree of freedom, 0, except 
in the special case, a = 0. For this special case, one has 
to distinguish two situations. If the bottom spins of a 
unit are antiparallel, there are two degrees of freedom 
when choosing the ground state orientation of the upper 
pair of spins: when counting ground-state coordinates, 
the loss of a coordinate that follows from the additional 
constraint, a — 0, is exactly balanced by the gain of an 
additional degree of freedom for the upper pair of spins. 
In the alternative situation, in which the first pair of 
spins is parallel, no freedom remains, and we obtain a 
lower-dimensional submanifold of the ground state. 

At this stage, we have fixed the orientation of spins 
in the lowest two layers of the lattice. Repeating the 
procedure in the subsequent layers, the value for a in each 
unit is determined by earlier choices, while one degree of 
freedom, <f>, remains for each unit of the system. 

Since any given ground state can be built up (or 
copied) layer by layer in this way the construction can be 
used to generate all possible ground states. By this con- 
struction we have demonstrated that the extensive part 
of the dimension of the ground-state manifold is equal to 
the number of units in the system. 



We show that the ground-state manifold is connected 
by demonstrating that any ground state can be continu- 
ously deformed into any other ground state without cost 
in energy. To do so, we choose a reference ground state 
and give an explicit construction by which the reference 
state can be reached from any ground state, without leav- 
ing the ground state manifold. This is done by consid- 
ering successive layers of the system, and bringing the 
spins in each layer into the orientation that they have in 
the reference state, using moves within the ground-state 
manifold which leave spins unchanged in the layers al- 
ready visited. The reference state is chosen to be one in 
which all spins on lattice sites equivalent under transla- 
tion have the same orientation, and in which each spin is 
antiparallel to the other spin belonging to the same layer 
and unit. 



o spins to be moved 
• spins held fixed 



step A 



step B 
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FIG. 5. When altering the orientation of the spins repre- 
sented by empty circles, only spins inside the wedge defined 
by the fixed spins (filled circles) have to be adjusted. For the 
pyrochlore lattice, a cone replaces the wedge. 

We make use of the following facts, which follow 
straightforwardly from consideration of the ground-state 
configurations illustrated in Fig. |]: 

(I) For a single unit of four spins, the manifold of 
ground states is connected: a and <f> (Fig. [|), which pro- 
vide a complete parametrisation of the internal degrees 
of freedom of the ground-state manifold, can be chosen 
independently from intervals of the real axis. Moreover, 
the orientation of any two chosen spins can be changed 
arbitrarily and continuously without leaving the ground- 
state manifold, as long as the other two spins are uncon- 
strained. 

(II) If two spins in a unit are antialigned, so will be the 
other two, whose common axis can then be rotated ar- 
bitrarily and continuously while keeping the orientation 
of the first pair fixed and remaining in the ground-state 
manifold. 

(III) It is possible to change the orientation of one spin 
in a unit arbitrarily and continuously while that of a sec- 
ond is held fixed, without leaving the ground-state man- 
ifold, as long as the other two spins are unconstrained. 
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This is a special case of (I). 

(IV) From II and III it follows that one can continu- 
ously change the orientation of a pair of spins belonging 
to a unit in the bottom layer, or a pair of antiparallel 
spins in a higher layer, and remain within the ground- 
state manifold of the whole system, whilst keeping fixed 
the spins that lie outside a wedge (or cone) for the square 
lattice with crossings (or the pyrochlore lattice), as de- 
picted in Fig. |[ 

As a consequence, we can again work through the lat- 
tice layer by layer. The procedure is as follows. 

(A) Align the spins in the bottom layer, unit by unit, 
to coincide with the reference state. This fixes the two 
spins of each unit in the bottom layer to be antiparallel. 
As we adjust the orientation of spins in the lower layer of 
each unit, those in the upper layer of the same unit can, 
by I, be reorientated to keep the unit always in a ground 
state. At the same time, by IV, spins in higher layers can 
also be reorientated to keep the system as a whole in a 
ground state. 

At this stage, spins in the bottom layer are in their 
reference state. 

(B) Spins in the second layer now form antiparallel 
pairs, since they belong to units which have antiparallel 
pairs in the lowest layer. The spins in the second layer 
can therefore, by II, be adjusted to coincide with those 
in the reference state. While this is done, by IV, spins 
in higher layers can be concurrently reorientated to keep 
the system within a ground state, as in step A. 

From the way we chose the reference state, it now fol- 
lows that all neighbouring spins in the second layer are 
pairwise antiparallel. Therefore, we can repeat step B for 
the third and all higher layers. Once we have done so, all 
spins in the system coincide with those in the reference 
state. Since one can go between two arbitrary ground 
states via the reference state, this completes our proof 
that the ground state manifold is connected. 

C. Pyrochlore antiferromagnets with general n > 2 

The arguments presented in the previous section for 
Heisenberg antiferromagnets generalise directly to n- 
component spins with n > 3. For the ground-state con- 
struction, the main difference is that spins 2 and 3 in 
Fig. ^ can now be rotated in n — 2 directions, so that 
the extensive part of D is N(n — 2), as expected from 
Maxwellian counting. 

The construction of ground states for the XY model 
is much simpler than for the Heisenberg model: as in 
the case of a Heisenberg antiferromagnet on the kagome 
lattice, for a generic ground state the number of con- 
straints is equal to the total number of degrees of free- 
dom. The construction of the ground states therefore 
involves fewer choices. It proceeds as follows. The ori- 
entations of spins in the bottom layer can be chosen ar- 
bitrarily. The ground-state configurations of spins in the 



next layer are then almost completely determined, since 
each tetrahedron has two pairs of antiparallel spins. The 
only freedom remaining is the discrete choice of which 
spin to place on which of the two sites of the unit in the 
next layer, unless the two spins of a unit in the lower 
layer happen to be antiparallel. In this case the orienta- 
tion of the pair (which has to be antiparallel) in the top 
layer be chosen freely. 

The proof of the connectedness of the ground-state 
manifold, presented above, follows essentially from the 
connectedness of the ground-state manifold of a single 
unit and from the fact that for any orientation of a pair 
of spins in a unit, the other pair can be chosen so that 
the total spin of the unit vanishes. This, along with the 
other steps, carries over to the case n > 2. 

D. Ground-state correlations of the Heisenberg 
pyrochlore antiferromagnet 

It is clear in our construction of ground states that 
spin correlations are not propagated efficiently. In this 
section, we show that, nonetheless, a few long-range cor- 
relations in high-symmetry directions are built into the 
ground states of the pyrochlore Heisenberg antiferromag- 
net. We discuss the signatures of these correlations in 
magnetic neutron diffraction. 

1. Correlations between planes 

From the ground-state condition, L Q = 0, it follows 
that the sum of all the spin vectors in two adjacent (100) 
planes is zero (this sum is also the sum of L Q over tetra- 
hedra making up the two planes). Therefore, adjacent 
planes are antiferromagnetically correlated. Since these 
correlations are long ranged, we expect sharp peaks in the 
neutron scattering cross section in the (200) directions. 

These peaks differ from Bragg peaks in two ways. 
First, their amplitude scales differently with sample size. 
Consider a sample in the form of a cube of side L, and 
let the total magnetisation of a (100) plane be M. In 
a typical ground state, M ~ L and the peak scattering 
amplitude varies as M ■ L ~ L 2 , in contrast to L 3 for a 
Bragg peak. Second, they are sharp in only one direc- 
tion in reciprocal space. Consider scattering at a point 
displaced from (200) by the vector (qn,q±). We argue 
that, at fixed qj_, the scattering amplitude as a function 
of <j 1 1 has a peak centered on gii = 0, of width Sqi\ cx |q_i_|. 
Contours of constant scattering intensity therefore have 
a distinctive bow-tie shape. To understand in detail the 
reason for this, it is necessary to examine the correlation 
in the magnetisation, M, of a region of a (100) plane 
with linear size £, and that of its equivalent, displaced 
by a distance z in the [100] direction. Let 5Ai z be the 
difference between these magnetisations, and measure z 
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in units of the plane spacing. The magnetisation differ- 
ence for adjacent planes, 5Mi, arises entirely from spins 
belonging to tetrahedra that are only partially included 
in the region. The number of such spins is proportional 
to the size of the boundary of the region - and hence to C. 
Since there are only weak correlations between individual 
spins on distances larger than the size of a tetrahedron, 
we obtain 5M\ oc \/~C. Increasing z, 5M. Z follows a ran- 
dom walk: 5M. Z oc \fz-\fL. Since M. ~ C, we obtain a 
correlation length £ oc £, and therefore Sq» oc |q_i_|. 

Similarly, there are long-range correlations in the (111) 
directions. The {111} planes are alternately kagome and 
triangular planes. Adjacent kagome planes contain the 
bases of adjacent tetrahedra, which, in the intervening 
triangular planes, share a common apex. In any ground- 
state, the total magnetic moments of all (100) kagome 
planes are equal, and also opposite to the total magnetic 
moment of all (100) triangular planes. 



2. Consequences for neutron scattering experiments 

In recent neutron scattering experiments on., a.. single 
crystal sample of CsNiCrF§ by Harris et aZ.,E£E30 the 
angular dependence of the neutron scattering cross sec- 
tion is studied. The correlation length is longest in the 
(100) direction, shortest in the (110) direction, and inter- 
mediate in the (111) direction. Since the presence of two 
species of magnetic ions in CsNiCrFg makes a detailed 
comparison with theprv-difhcult, Harris et al. also report 
Monte Carlo studiesEMZl of the angular neutron scatter- 
ing cross section in the [hhl] plane, shown in Fig. [j[ which 
is taken from Ref pq . 

The scattering is broad in most directions except the 
[100] and [111]- directions, where narrow necks appear 
at low temperature. The scattering near these necks has 
the appearance of a bow-tie, as described above. In addi- 
tion to these bow-ties, which are narrow in the direction 
parallel to the wavevector transfer, there are subsidiary 
bow ties, narrow in a perpendicular direction. Their ori- 
gin can also be explained using arguments of the kind 
described above. 




FIG. 6. Static neutron scattering cross section in the 
[hhl] plane for a Heisenberg pyrochlore antiferromagnet, ob- 
tained numerically by quenching the system into a ground 
state. Horizontal and vertical coordinates are I and h respec- 
tively. Light shaded areas represent high intensity, and black 
represents zero intensity. 



E. Local zero modes 

We next discuss the existence of degrees of freedom in 
the ground state which involve only spins in a finite re- 
gion in the bulk of the system. These we call local zero 
modes. They are of interest because unhindered rotation 
of finite numbers of spins is likely to be particularly ef- 
fective at destroying correlations, both in space and in 
time. 

In the kagome magnet, the nature of the zero modes 
was established in Refs. ^ and The requirement 
that L = leads the three spins in each triangle to be 
coplanar and at relative angles of 120°. If all the spins on 
the lattice are coplanar, there are three spin orientations, 
(A, B,C), each of which occurs in each triangle once. 
A zero mode, called a weathervane defect,E3 arises as 
follows. In a region enclosed by a line of spins of one type 
(say, C) , a line of spins (alternating between type A and 
B) can be rotated about the spin direction of C at no cost 
in energy. Starting from a particular state, gray, ground 
state can be constructed using these zero modes.EaH 

There is a closely related way of describing the zero 
modes for the XY magnet on the pyrochlore lattice. 
They are again associated with closed loops, this time 
of antiferromagnetically oriented spins. Our construc- 
tion of a ground state in section III amounts to finding a 



set of lines through the lattice sites with adjacent sites on 
a line occupied by antiparallel nearest neighbour spins, 
and each site belonging to exactly one line. Each line 
can be labelled with an angle giving the orientation of its 
spins, and a zero mode involves changing one such angle. 
VillainEij used a description of this kind to generate a 
subset of the ground states of the Heisenberg pyrochlore 
antiferromagnet. By contrast, for the XY model, this 
approach generates all the ground states. 

For the Heisenberg antiferromagnet on the pyrochlore 
lattice, the nature of the zero modes is much more com- 
plicated, reflecting the larger freedom in the ground-state 
manifold. We have not been able to find a simple descrip- 
tion of a generic zero mode. We nevertheless believe that 
in a sufficiently large region of a generic ground state, 
there are local zero modes. Our argument rests on count- 
ing degrees of freedom: the number of degrees of freedom 
a region contributes to the ground state is proportional 
to its volume, while fixing surrounding spins imposes a 
number of constraints proportional to its surface. For a 
large enough volume, the number of degrees of freedom 
exceedes the number of such constraints. The existence 
of such local modes in all ground states is however not 
guaranteed. For instance, for a state in which all spins in 
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a given (001) plane are parallel, and antiparallel to the 
spins i n neig hbouring planes, the sum rules discussed in 
section [ill D| preclude the existence of local zero modes. 



IV. GROUND-STATE SELECTION AT LOW 
TEMPERATURES: ANALYTICAL RESULTS 

In this section, we examine the circumstances under 
which thermal fluctuations induce order in geometrically 
frustrated antiferromagnets. This phenomenon - known 
as order by disorder - has been discussed in great detail 
for the Heisenberg antiferromagnet on the kagome lat- 
tice, where th^pap^ fluctuations induce coplanar ordering 
of the spins .cJc3'E2l There is evidence from past simula- 
tions that order by disorder is not a universal occurrence 
in geometrically frustrated systems - both Heisenberg 
spins on the pyrochlore latticeEjO and four-component 
spins on the kagome latticeEil apparently remain disor- 
dered at low temperature - but the systematics have not 
previously been studied. 

The experimental situation for pyrochlore magnets is 
rather complicated. A few compounds develop long- 
range order at low temperatures ,E3 whereas others un- 
dergo a spin-glass transition.Eij It is unclear what the im- 
portance is of various features of real systems. Neverthe- 
less, we restrict our attention in the following to the theo- 
retically idealised problem of a classical antiferromagnet, 
without anisotropy, disorder, further-neighbour or dipo- 
lar interactions. 

This material is arranged as follows. First, we consider 
the analytically accessible problem of four Heisenberg or 
XY spins on a single tetrahedron. Next, we investigate 
the general case of a lattice built from groups of q spins, 
each with n components, and ask whether thermal fluc- 
tuations restrict the spins to an ordered (e.g. collincar 
or coplanar) configuration. Then, in Sect. M, we present 
the results of numerical simulations, which test the con- 
clusions reached from our analytical arguments. 



A. The single tetrahedron 

We first study a problem simple enough to allow ex- 
plicit evaluation of some of the quantities of interest: an- 
tiferromagnetically coupled spins occupying the corners 
of an isolated tetrahedron. 

There are eight degrees of freedom associated with 
four Heisenberg spins, and for the system to be in a 
ground state, three constraints must be satisfied, since 
L = 0. The ground state manifold therefore has five 
dimensions: of these, three arise from global rotations, 
while the remaining two can be parameterised as dis- 
cussed in Sect. Ill A. The energy cost for fluctuations 



special ground states in which spins are collinear, energy 
varies quadratically with displacement from the ground 
state manifold only in two directions, and quartically in 
the third, indicated schematically in Fig. [?]. The collinear 
states are the obvious candidates for selection by thermal 
fluctuations. 




FIG. 7. Soft fluctuations around a collinear state 

To study such selection, we have calculated the prob- 
ability distribution, W(8), for the angle, 6, between a 
pair of Heisenberg spins, integrating over all orienta- 
tions of the four spins wi th a Boltzmann distribution and 
the Hamitonian of Eq. 2.1. Two factors contribute to 



from this ground-state manifold in the remaining three 
directions in configuration space is, for a generic ground 
state, quadratic in displacement. By contrast, for the 



W(6)d6: the measure, sin(6)d6, and a statistical weight. 
The low-temperature limit of the latter is [2 cos((9/2)] -1 , 
the divergence as 9 — > 7r reflecting the lower free en- 
ergy attached to fluctuations around the collinear state. 
Combining both factors, W{6) = sin(0 /2)d0 /2: config- 
urations which are nearly collinear have higher weight 
than others in this distribution, but the entire ground- 
state manifold is accessible even in the low-temperature 
limit, and there is no fluctuation-induced ground-state 
selection. 

To illustrate the alternative, consider the same prob- 
lem for XY spins. In this case, the measure contributes 
simply dd to the distribution, W(0)dd, while the sta- 
tistical weight at temperature T <C J is proportional 
to Isin^)]- 1 for sin 2 (6») > T/J, and to (J/T) 1 / 2 for 
sin 2 (#) <C T/J. As a result, the weight in the limit 
T — > is overwhelmingly concentrated near collinear 
spin arrangements (9 = and 7r), reflecting selection 
of these states by thermal fluctuations. Order by dis- 
order is just such a concentration of statistical weight on 
a submanifold of ground states. Note that it can occur 
in a finite system (in this case, a system of four spins), 
and is quite different from the order that appears in a 
symmetry-breaking phase transition, which is restricted 
to the thermodynamic limit. 

It is straightforward to demonstrate these effects in 
Monte Carlo simulations. In Fig. [§, we pl ot the collincar- 
ity parameter, P(l) (defined in section V A , Eq. |5.l| ), 
as a function of Monte Carlo time, for a simulation of 
four Heisenberg spins arranged in a single tetrahedron 
at the temperature T = 2.5 x 10 _5 J. For the current 
purposes, it is sufficient to note that the collinearity pa- 
rameter takes on values between —1/3 (for a state with all 
spins at relative angles of 70.5° or 109.5°) and +1 (when 
all spins are collinear). We see from Fig. |^ that the sys- 
tem explores all ground states, attaining values of the 
collinearity within 2 x 10~ 4 of the extremal ones, and is 
not trapped near a collinear state. The average of the 
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collinearity parameter, 0.193 ± 0.02, is distinct from 0, 
its value in the high-temperature limit, and close to the 
exact low-temperature value, 1/5, obtained from the ex- 
pression for W{&) given above. 




FIG. 8. Evolution of the collinearity parameter with 
Monte-Carlo time for a single tetrahedron at T = 2.5 x 10~ 5 J 
for 2 x 10 Monte-Carlo steps per spin. 



B. The general problem: groups of q spins with n 
components 

We now examine whether thermal fluctuations select 
particular ground states for the general class of sys- 
tem introduced in Sect. |l[ in which a lattice is built 
from N corner-sharing units of q spins, .each having 
n components. We have argued elsewhere^ that low- 
temperature behaviour is characterised by a probability 
distribution over the ground-state manifold, defined in 
the limit T — > 0. Let x be coordinates on the ground- 
state manifold; at each point x, one can introduce local 
coordinates, y, spanning the remaining directions in con- 
figuration space. Generically, the energy of the system 
relative to its ground-state value will have a Taylor ex- 
pansion with the leading term 



resulting in a ground-state probability density 



Z(x) oc / d{yi} 



-PH 2 



Il^BT/eKx)) 



1/2 



(4.1) 



(4.2) 



In principle, order might arise in cither of two ways. 
First, it can happen that certain, special ground states 
have soft fluctuations, so that at some point, Xo, on the 
ground state manifold (or, more generally, on some sub- 
space), some of the ej(xo) vanish. Then Z(x) will di- 
verge as x approaches xo. If any such divergences are 
non -int egrable, one should keep higher order terms from 
Eq. 4.1 when calculating Z(x). The result of doing so will 
be, in the limit T — ► 0, a distribution concentrated exclu- 
sively on the subset of ground states for which Z{x) is di- 
vergent: these are the configurations selected by thermal 
fluctuations. It is this mechanism for fluctuation-induced 
order that we study. There is, however, also a second pos- 
sibility, which we do not persue here: it might happen 
that the probability density, Z(x), is spread smoothly 



over the ground-state manifold, but that there never- 
theless exist correlation functions which, when averaged 
with this weight, are long-ranged. 

To decide whether ground states with soft modes are 
selected, it is necessary to know the number, M, of e; that 
vanish, and the dimension, S, of the subspace on which 
this happens. Close to this subspace, we separate x = 
(u, v) into an S'-dimensional component u, lying within 
the subspace, and a [D — S')-dimensional component v, 
locally orthogonal to it, with magnitude v. We expect at 
small v the behaviour e;(x) oc v 2 for M of the e/'s. Hence, 
Z{x) diverges as v~ M for small v, and the subspace is 
selectedEil as T — > if the integral 



J Z(u, v)dv oc J v 



D-S-M-l 



dv 



(4.3) 



is divergent at small v. 
n 



Disorder 




3 4 5 6 <! 

FIG. 9. The occurence of order by disorder for 

n-component spins arranged in corner-sharing units, each con- 
sisting of q spins. Ordered (marginal) models are denoted by 
circles (crosses). 

We therefore need to consider candidate ordering pat- 
terns, and determine the sign of D — S — M in each 
case. It seems in general that the preferred ground 
states are ones in which spins are collinear or coplanar, 
because these have the largest number of soft modes. 
Collinear spin order is possible on lattices built from 
units containing an even number of sites, q: in practice, 
those constructed from tetrahedra. Such order results 
in one soft mode per unit, as illustrated in Fig 0. The 
number of soft modes is therefore M — N, and (since 
S = n — 1) we expect order only if D < N. Estimating 
D as F — K = N[q(n — 1) — 2n]/2, we predict order if 
n< (q + 2)/(q - 2), and disorder if n > (q + 2)/(q - 2). 
Thus, for pyrochlore antiferromagnets, two-component 
spins order, and four-component spins do not. The ap- 
proach reaches no conclusion in the marginal case of 
three-component spins, but simulations (Refs. |3^,^8| and 
as described below) indicate that Heisenberg spins do not 
order. On lattices made from corner-sharing triangles, 
such as the kagome lattice, there are no collinear ground 
states; instead, coplanar ordeip4jnay occur. Such order re- 
sults in M = N/2 soft modesjHJ using again the estimate 
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F — K = N[q(n — 1) — 2n]/2 for D, we predict order in 
this case if n < 4 and disorder if n > 4. Simulations of 
kagome antiferromagnets demonstrate that there is in- 
deed coplanar order for n-^F 3 £3 and that the marginal 
case, n = 4, is disordered. Eil 

Summarising, the only cases in which there is order by 
disorder are q = 4, n — 2 (the pyrochlore model) 
and q = 3, n = 3 (the Heisenberg kagome model). In 
both instances, there is a low entropic cost to enter the 
ordered state (Dm = 0) and a high entropic gain from 
soft fluctuations because in the ordered state the con- 
straints, L a = 0, are not independent. These conclusions 
are depicted in Fig. |j. 



V. GROUND-STATE SELECTION AT LOW 
TEMPERATURES: NUMERICAL RESULTS 

In the following two subsections, we present the re- 
sults of the Monte Carlo simulations of XY and Heisen- 
berg pyrochlore antiferromagnets. One aim is to test our 
prediction of collinear ordering for XY spins. We also 
consider the Heisenberg model in detail, to show that 
order by disorder is indeed absent in this case. Our stud- 
ies of the Heisenberg antiferromagnet are a continuation 
of Reimers' pioneering simulationsE3 and Zinkin's sub- 
sequent work.E3 Our conclusions are in agreement with 
these authors, in particular with the earlier - albeit ten- 
tative - ideas of Zinkin, but our results are more ex- 
tensive. Reimers work concentrated on the temperature 
range T > 0.05J: many of the observations described in 
the following are very hard to discern or absent in this 
regime. 

Our simulations were carried out on systems of sizes 
ranging from one unit cell (N = 2 tetrahedra, N s — 
4 spins) to 17 3 unit cells (N = 9826, N s = 19652). 
As pointed out in section IV, small systems display 
large fluctuations, and therefore require very long sim- 
ulation runs. For N s — 4, the longest simulation was 
2 x 10 8 Monte Carlo steps per spin at T = 5 x 1CT 5 J. 
For the largest system, however, only 1.5 x 10 6 Monte 



Carlo steps per spin were necessary even at the lowest 
temperature. 



A. Correlation functions 

In this subsection, we consider two-spin correlations 
and also a correlation function which quantifies directly 
the collinearity of the spin system. In the next subsec- 
tion, we discuss the heat capacity, which is an indirect 
probe of the state of the system, but in some ways more 
conclusive, since it is sensitive to the presence of soft fluc- 
tuations irrespective of the type of ordering with which 
they are associated. 

First, we demonstrate that the Heisenberg model does 
not have Neel order, even at low temperature. The corre- 



lation function Q(r) = (S(0) ■ S(r)) is shown in Fig. |T^: 
correlations are very small beyond the second neighbour 
distance. Second, to measure the collinearity of spins, we 
evaluate the correlation function (for n-component spins) 



P(r) = 



1 



((S(O).S(r))' 



(5.1) 



which is constructed to have the values P = at infi- 
nite temperature and P = 1 in a collinear state. P(r) 
is shown in Fig. [l^, with r in units of nearest-neighbour 
distances. The correlations for Heisenberg spins again 
have a range of only two nearest neigbour distances: 
there is no fluctuation-induced order. Equally, the pre- 
dicted collinear order for XY spins is confirmed: there is 
long-range order in P(r) at this temperature. Note that, 
despite the very low temperature, the order parameter, 
P(r — ► oo ) ~ 0.86, is appreciably less than its maximum 
possible value of 1. We expect on general grounds that 
such nematic order should be established via a first-order 
phase transition, but have not attmepted to check this in 
detail in our simulations. 
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FIG. 10. Correlation functions for the Heisenberg and XY 
antiferromagnets at a temperature of T — 5 x 10~ 4 J. The 
two-spin correlation function, Q(r), (dot-dashed line) and the 
collinearity correlation function, P(r), (solid line) for a system 
of 2048 Heisenberg spins; and P(r) for a system of 864 XY 
spins (dashed line). 

The temperature dependence of collinearity for neigh- 
bouring spins is shown in Fig. |ll|. Neighbouring Heisen- 
berg spins have a limiting low-temperature value, P(l) — 
0.2, which is non-zero because the correlation length, 
though small, is itself finite. By contrast, XY spins 
become perfectly collinear in the low-temperature limit. 
The low-temperature variation of [1 — P(l)], the devia- 
tion of collinearity from its maximal value is character- 
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istic of fluctuat ion- in duced order. Specifically, we ex- 
pect 1 — P oc yT/J at low temperatures, because quar- 
tic modes give rise to the dominant fluctuations at low 
temperatures. These modes, with coordinates 77, charac- 
terised schematically in Fig. ^, have 77 oc T 1 / 4 by equipar- 
tition. Since (Si • S 2 ) 2 - (1 - rj 2 /2) 2 ~ 1 - rr 2 , we obtain 
1 - P(l) oc T 1 / 2 . We show in Fig. |l| that P(l) does 
indeed behave in the expected way. 

We have checked the dependence of our results on 
length of simulation run and system size. To test whether 
the system is properly equilibrated during our Monte 
Carlo runs, we investigate the dependence of data on 
initial conditions, comparing results from random and 
collinear initial states. For Heisenberg spins, our sim- 
ulations are long enough that neither of the correlation 
functions studied retains memory of the initial state. For 
XY spins, we are able to equilibrate P(r) (see Fig. 11), 
but not Q(r): collinear order presumably hinders relax- 
ation of two-spin correlations. To test for finite-size ef- 
fects, we carry out simulations on systems ranging in size 
from N s — 4 to N s — 19652 spins. Only for small sys- 
tems of XT-spins are marked finite size effects observed, 
as shown in Fig. uji. 
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FIG. 12. P(r = 1) versus (T/J) 1/2 for XY spins. 
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FIG. 11. The temperature dependence of P(r = 1) for XY 
and Heisenberg spins. Most error bars are smaller than the 
symbols. Simulations started from random initial spin config- 
urations, except those for the data points marked with open 
circles, which started from collinear spin configurations. 
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FIG. 13. System size dependence of the collinearity param- 
eter P(r — 1) for XY and Heisenberg spins at T = 5 x 10~ J 
and T = 5 x 10" 4 J. 



B. Specific heat 



An unbiased way to search for soft fluctuations and - 
by implication-^- fluctuation-induced order is to measure 
heat capacity.EfJ At low temperature, one expects to be 
able to describe fluctuations of the system from a ground 
state in terms of canonical coordinates which are almost 
independent of each other. From the classical equiparti- 
tion theorem, a canonical coordinate, p, which appears in 
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the Hamiltonian as (p/po) 2r , contributes A;s/(2r) to the 
heat capacity. Hence, each quadratic mode contributes 
ks/2, and each quartic mode fcs/4, whereas zero modes 
do not contribute at all. Determining the heat capac- 
ity therefore allows a determination of the number of 
quadratic and quartic modes present. 

If, as predicted, thermal fluctuations select a collincar 
state for the XY model, then there is one quartic and one 
quadratic mode per tetrahedron. This results in a heat 
capacity per spin, C, of 3/8ks- In the absence of order, 
all modes are quadratic, and C — ks/2. For the Heisen- 
berg antiferromagnet, there are four degrees of freedom 
per tetrahedron, one of which is a zero mode. If there 
is no order, we expect C = 3/4fce; if a collinear state is 
selected by thermal fluctuations, C — 5/8fcg. In finite- 
sized systems, the heat capacity per spin is reduced. For 
our choice of periodic boundary conditions, this mani- 
fests itself in the correction C(N) = ((N - 1)/N) C(oo). 

As shown in Fig. |lj, we find in the limit N — > oo, 
that C = 0.376 ± 0.002 for XY spins, and that C 
!i.7 17 ± 0.002 for Heisenberg spins. This is consistent 
with the presence of order for XY spins, with one quar- 
tic mode per tetrahedron, as expected. For Heisenberg 
spins, we obtain an upper limit of 0.04 quartic modes per 
tetrahedron. 



0.75 




FIG. 14. The specific heat for different system sizes for the 
XY and the Heisenberg antiferromagnet on the pyrochlore 
lattice. The lines are a guide to the eye. 



VI. THE DYNAMICS OF THE PYROCHLORE 
ANTIFERROMAGNET 

We now turn to time-dependent correlation functions, 
and ask how the system explores the vicinity of its 
ground-state manifold at low temperature. We study the 



spin autocorrelation function with precessional dynam- 
ics, both analytically and by numerical integration of the 
equations of motion. Two of the facts established in the 
previous sections have important implications for spin 
dynamics. First, we have shown that the ground-state 
manifold is connected, which means that the magnet docs 
not get trapped in a particular state at low temperatures 
by internal energy barriers. Second, we have shown that 
there is no entropic selection of special ground states at 
low temperatures, which suggests that the dynamics is 
unlikely to be hindered by free energy barriers. We in- 
deed find that spin correlations relax relatively rapidly 
even at low temperatures, having a time-scale that di- 
verges only as T _1 , and not, for example, according to 
an Arrhenius law. 



A. Derivation of an effective spin dynamics 

Consider, in the first instance, the dynamics linearised 
around a ground state. In this approximation, a sys- 
tem of N tetrahedra (hence with 2N spins and a 47V- 
dimensional phase-space) will have 2N normal modes. 
Some of these modes will be conventional, finite fre- 
quency spin- waves, but a fraction will have zero fre- 
quency, because there is no restoring force for displace- 
ments in phase space that take the system from one 
ground state to another. Beyond the harmonic approxi- 
mation, non-linear terms in the full equations of motion 
will have various consequences: the conventional, finite- 
frequency modes will acquire a finite lifetime; and cou- 
pling between these modes and the ground-state coordi- 
nates will drive the system around its ground-state man- 
ifold. We find that there are three distinct time-scales 
at low temperature. The period of the highest frequency 
spin-waves, 0(K/J), sets the shortest scale; their lifetime, 
t s ~ ^/[J/cbT] 1 / 2 , provides an intermediate scale; while 
the longest scale is the decay time of the autocorrela- 
tion function, r ~ h/hsT. This separation of time-scales 
greatly simplifies the problem. 

Our starting point is the equation of motion, 



f = S,xH,( t ) 



-JS, x (L Q + l>p), 



(6.1) 



where we have set ft = 1. Hj(t) is the exchange field 
acting at site i, which can be expressed in terms of L Q 
and Lg, the total spins of the two tetrahedra to which 
Si belongs. Summing over the sites of a tetrahedron, the 
time-dependence of L a is 

"' L " -J^S^xL^, (6.2) 





dt 



in which the notation S Q ^ has been introduced for the 
spin common to the tetrahedra a and (3. 

The right side of Eq. 3.2 implicitly defines a 37V x 37V 
matrix, M, acting on a vector constructed from the com- 
ponents of the L a 's. This matrix, being real and anti- 
symmetric, has eigenvalues which are purely imaginary 
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and occur in pairs, ±iui, related by complex conjuga- 
tion. For a ground-state spin configuration, the magni- 
tudes of these eigenvalues are the frequencies of 3N/2 
normal modes, while the real and imaginary parts of the 
associated eigenvectors are canonically conjugate coordi- 
nates for the modes. The remaining N directions in phase 
space are spanned by coordinates having L Q = for all a, 
and therefore lie within the ground-state manifold. The 
matrix M is well-defined and has purely imaginary eigen- 
values for any spin configuration; for a low-temperature 
spin configuration, the eigenvalue magnitudes are pre- 
sumably good approximations to the normal mode fre- 
quencies in a nearby ground state. We display in Fig. [l5| 
the density of states, p(tv), on a linear scale, obtained 
by diagonalising M for low-temperature pyrochlore spin 
configurations generated in a Monte Carlo simulation. It 
is noteworthy that p(w) appears to be finite at lo = 0: 
p(uj) neither includes a divergent contribution, propor- 
tional to S(ui), nor does it vanish as u> — > 0. 
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FIG. 15. The density of states, p(ui) for a system of 2048 
Heisenberg spins at low temperature. 

The fact that p(uS) does not contain a delta function 
at io — gives information on how canonically conju- 
gate pairs of coordinates appear in the linearised dynam- 
ics. Quite generally, the Hamiltonian in the harmonic 
approximation can be reduced to the form 



H 



2iV 

E 

i=i 



(aipf + faq?), 



(6.3) 



where pi and qi are a canonically conjugate pair of co- 
ordinates. We know from the ground-state construction 
described earlier that 1/4 of these coordinates belong to 
the ground-state manifold, and therefore that 1/4 of the 
AN numbers {a;, f3i} are zero. Oscillations of the coor- 
dinates pi, qi have zero frequency if ai ■ j3i = 0, and so 
the fraction of zero- frequency spin- wave modes might, in 
principle, range from 1/4, if for each of these modes both 
cti = and j3i = 0, to 1/2, if for each of these modes only 



one of ai and j3i is zero. The fact that p(u>) does not con- 
tain a delta function at u) = implies that all (except for 
a fraction vanishing in the thermodynamic limit) of the 
3./V/2 modes derived from M have non-zero frequency, 
and therefore that only the remaining 1 /4 of modes have 
zero frequency. Hence, for every I, either both ai and (3i 
are zero or neither is zero: coordinates in the ground- 
state manifold all appear in canonically conjugate pairs. 

The non-zero density of states apparent at small 
frequency in Fig. |l5| is in striking contrast to the 
behaviour, p(u>) oc uj 2 , that occurs in both Neel- 
ordered antiferromagnets (including ordered states of 
the pyrochlore antiferromagnetEj) and conventional 
spin-glasses. The arguments used by Halperin and 
Saslowo and by GinzburglEI to predict propagating long- 
wavelength modes in spin-glasses depend on the ground 
state having a stiffness. This stiffness appears to be miss- 
ing in ground states of the pyrochlore antiferromagnet, 
because of the many zero modes. 

The exchang e fi eld, Hj(t), appearing in the equation 
of motion, Eq. 6.1, can be written as a superposition of 



contributions arising from the finite-frequency modes, in 
terms of vectors, u;(i), determined by M, and ampli- 
tudes, Ai 7 determined by the initial conditions: 



Hi(t) 



l 



Am{i)e Wlt + c.c 



(6.4) 



In the harmonic approximation, the amplitudes, Ai, are 
time-independent, but in the full dynamics their magni- 
tude and phase will change on a time-scale which defines 
the spin- wave lifetime, t s . We postpone detailed discus- 
sion of the temperature dependence of t s until the end 
of this section, but note that, at low temperature, r s is 
large compared to the typical spinwave period, J -1 . 

If the equation of motion is integrated over time- 
intervals longer than t s , contributions from modes with 
frequencies uji 3> t s _1 average to zero, while those from 
modes with u>i < t^ 1 fluctuate randomly, according to 
the time-dependence of Ai(t). Hence Hj(t) has mean 
value zero, and has fluctuations which are characterised 
most importantly by their low- frequency spectral density, 



In terms of the amplitudes, Ai(t), 

r = E [(M0)At(t))e iUlt \Mi)\ 2 dt. 



(6.5) 



(6.6) 



Since, for large r s , only the low- frequency modes con- 
tribute to r, and since, from equipartition, (|H;(i)| 2 ) ~ 
J 2 (|L a | 2 ) ~ Jk B T, we have T oc Jk B Tp(0+). Note that, 
as p(ui) oc J -1 , r is independent of J. 

We now pro ceed to calculate long-time spin correla- 
tions from Eq. 6.1, by treating Hi(t) as Gaussian white 
noise with the correlator 
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{U i (t)-H i (t'))=2r6(t-t'). 



(6.7) 



To do so involves several assumptions. The most impor- 
tant of these is that the spinwavc lifetime, t s , which sets 



the width in time of the delta function in Eq 6/7 , is small 
compared to the decay time of the spin autocorrelation 
function. We show below that this is asymptotically ex- 
act as T/J — > 0. Further, in taking T to be a constant, 
rather than a functional of the instantaneous spin con- 
figuration, we implicitly neglect variations with spin con- 
figuration in the local density of states, J2i l u « W| 2 ^( w 0- 
Such fluctuations certainly exist, but seem from our nu- 
merical studies only to have a small effect on the form of 
the spin autocorrelation function. Solving the Langevin 
equation that results from treating Hi(t) in Eq. |6.l| as 
white noise, we obtain (Sj(0) • S-j(t)) = e _r * and hence 
(reinstating K) 



(S, ; (0) • Si(t)> = exp(-ck B Tt/h) , 



(6.8) 



where c is a dimensionless constant of order unity. There- 
fore, the autocorrelation time r = = h/(cksT). We 
emphasise again that, at T <C J//cb, it is T alone, and 
not J, which sets the scale for long-time dynamics. 

To complete this discussion, it is necessary to estimate 
the spinwavc lifetime, t s . There are two physical pro- 
cesses that contribute to r s . One, common to all an- 
tiferromagnets, is the anharmonic interaction between 
different finite-frequency modes, which here results in a 
lifetime varying as T~ 1 for small T. It is, however, over- 
whelmed by a second process, specific to systems with 
many ground-state degrees of freedom, in which finite- 
frequency modes are mixed by the motion of the system 
between different ground states. More formally, on time 
scales 3> ti/J, the matrix M is time-dependent. The 
linearised equations of motion, with time-dependent M, 
define an autonomous dynamical problem in which the 
instantaneous normal mode amplitudes, Ai(t), are time- 
dependent. The time-dependence of the matrix elements 
of M mixes amplitude, initially concentrated in a sin- 
gle mode, I, over all modes lying within a window of 
frequencies around u)%. From time-dependent perturba- 
tion theory, a fractional change, /, in matrix elements 
spreads amplitude over a frequency window whose width, 
Slu, forms a fraction / of the entire spinwave spectrum, 
so that Su> ~ fJ/h. And from our results for the spin 
autocorrelation function, the fractional change in matrix 
elements during a time-interval r s is / = k-QTT s /h. The 
spinwave lifetime is the time at which the frequency win- 
dow resulting from this / has width 5u> ~ 1/t s , and so 
we obtain r s ~ ^/[J/cbT] 1 / 2 . As required for the con- 
sistency of our arguments, at low temperatures this is 
indeed a much shorter timescale than that for the decay 
of the spin autocorrelation function. 

As the temperature is raised towards T ~ J, this sep- 
aration of timescales breaks down. The precession on 
the previously shortest timescale then becomes visible, 
and the autocorrelation decays initially as 1 — A(t) cx t 2 



rather than 1 — A(t) cx t. This is indeed observed in 
our numerical simulations described in the next section 

(Fig. hd. 



B. Molecular dynamics simulations 

In this subsection, we present results obtained from 
numerical simulations of the dynamics of the Heisenberg 
pyrochlore antiferromagnct. In these simulations we eval- 
uate the autocorrelation function 



A(t) = (S,;(0) • Si(t)) . 



(6.9) 



We find that A(t) decays exponentially in time. The sim- 
ulations confirm the predictions of the preceding section, 
namely that the timescale for the dynamics, r, varies as 
T -1 . Finally, we do not discover any sign of spin freezing 
even at temperatures as low as T = 5 x 10~ 4 J. 

We generate uncorrelated, thermalised initial configu- 
rations by Monte Carlo simulation, from which the equa- 
tion of motion, Eq. |6.1| , is integrated using a fourth- 
order Runge-Kutta algorithm. Related calculations for 
the kagome Heisenberg antiferromagnet have been de- 
scribed previously by Keren. o 

Some details of our procedure are as follows. We 
choose the integration time-step so that energy is con- 
served to at least one part in 10 8 . The temperature 
range of the simulations covers three orders of magni- 
tude, the lowest temperature being T = 5 x 10~ 4 J. The 
system sizes studied range from 32 spins to 2048 spins. 
There are marked finite size effects in the smaller sys- 
tems, which we believe result from all spins precessing 
together about the total magnetisation, M tot , of the sys- 
tem. Since |Mt ot | 2 ~ N, the precession rate varies as 
N^ 1 / 2 , and decreases rather slowly with increasing sys- 
tem size. For the results presented, we hasten conver- 
gence to the thermodynamic limit by adding the term 
JM 2 ot to the Hamiltonian, which constrains the total 
magnetisation to be independent of system size and near 
zero. As a result, values of the decay time, r, coincide 
for systems with 500 and 2048 spins. 



1. The functional form of A(t) 



From the analytic calculation presented in subsection 



VIA , we expect A(t) to depend on time and tempera- 
ture only through the combined variable, Tt. We show 
in Fig. [l6| A(t) as a function of this scaling variable, at 
various temperatures and over one and a half decay times, 
t, for a system of 2048 spins. The collapse of the data 
onto a single curve, at all except the highest temperatures 
(T/J > 0.1), is striking evidence in support of our an- 
alytic results. To demonstrate the accuracy with which 
the decay of A(t) at low temperature is exponential, and 
to indicate the magnitude of finite-size effects in our re- 
sults, we show in Fig [l?] data for A(t) on a logarithmic 
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scale, for runs starting from different initial configura- 
tions generated at the same temperature. 
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FIG. 16. The autocorrelation function as a function of the 
rescaled time, Tt. 
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FIG. 18. The decay time as a function of tempeature for 
three different system sizes. 
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FIG. 17. The decay of A(t) on a logarithmic scale for five 
different runs at the same temperature, T/J — 6 x 10~ 3 , in a 
system of 2048 spins. 

To examine quantitatively the temperature depen- 
dence of the decay time, r, we fit data at each tempera- 
ture to an exponential, exp (— t/r). The resulting values 
for r are displayed in Fig. [l|. In order to extract the 
temperature-dependence of r, we fit it to the power law 
t = AT^^SA Excluding temperatures T/J > 0.15, we 
obtain C = 0.998 ± 0.012 and A = 0.53 ± 0.04. This re- 
sult agrees with and confirms our prediction that £ = 1. 



C. Inelastic neutron scattering 

Inelastic neutron scattering provides the most detailed 
probe of dynamical correlations. We expect from the re- 
sults described above that diffuse inelastic scattering in 
the temperature range Tp < T -C |0cw| should have a 
Lorentzian lineshape in energy, with a width, F, varying 
as r = cksT, where c is of order unity. Although we have 
not explicitly examined the dependence of dynamic corre- 
lations on wavevector, k, it seems likely that main contri- 
bution to the inelastic linewidth in the low-temperature 
limit should be roughly wavevector independent. One 
reason for thinking this is that the decay of the autocor- 
relation function probably arises because of rotation of 
relatively small clusters of spins - the local ground-state 



degrees of freedom identified in III E . In consequence, we 
expect dynamical correlations to be short-ranged in space 
and broad in wavevector. In addition, conservation laws 
which might result in significantly different behaviour, 
for example T oc k 2 for small k from conservation of spin 
density, do not appear to be in operation: since the mag- 
netisations of individual tetrahedra are identically zero in 
classical ground states, instantaneous magnetisation fluc- 
tuations can decay without spreading to large distances. 
Thermally induced fluctuations in the magnetisations of 
tetrahedra may result in an additional, diffusive compo- 
nent to spin correlations, with an amplitude that vanishes 
in the low-temperature limit. 

Inelastic neutron scattering from the pyrochlore anti- 
ferromagnet|-C , siVzCri ; e, which has OcwH — 70K and 
Tp — 2.2K,E3 is reported in Ref. [l6, in which the en- 
ergy dependence of scattering is fitted by a Lorentzian. 
The linewidth decreases as temperature is decreased be- 
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low 70 K, but the data seem insufficiently precise to test 
whether T is linear in T. Similar experiments on SCGO, 
in which |0ow| is around 500 K and Tp ~ 3.5K, yield a 
Lorentzian inelastic lineshape with a temperature depen- 
dence of the width which is encouragingly close to linear, 
and of oi'|dp£, fc^T, over the temperature range from 30A" 
to 290A'JlI2 



VII. CONCLUDING REMARKS 

We have presented a detailed theoretical analysis of 
the low-temperature properties of a class of geometri- 
cally frustrated classical antiferromagnets, with particu- 
lar emphasis on the pyrochlore Heisenberg antiferromag- 
net, which has a macroscopically degenerate ground-state 
manifold. Given the connectedness of this ground-state 
manifold and the absence of appreciable free energy bar- 
riers, there seems to be no mechanism for localising the 
system in a particular region of the ground-state mani- 
fold, and it is therefore unlikely that the spin glass tran- 
sition observed in most experiments is a feature of the 
disorder-free classical isotropic Heisenberg pyrochlore an- 
tiferromagnet. Rather, we find correlation functions to 
be short-ranged in both space and time, and conclude 
that the spins continue to fluctuate strongly down to the 
lowest temperatures. 

We have analysed for the first time the low-energy 
dynamics of these geometrically frustrated antiferromag- 
nets. Our discussion does not depend on details of the 
pyrochlore lattice structure. In fact, we expect it also 
to apply to the Heisenberg model defined on the SCGO- 
lattice of Fig. || since that model has Dm /N > and 
does not - following the arguments presented in Ref. |5l] 
and Sect. [iy| - display order by disorder. These proper- 
ties are in striking contrast to those of the kagome Heisen- 
berg antiferromagnet, and SCGO is much more similar 
to a pyrochlore magnet than to the kagome system. 

We expect our results to be robust against-jthe intro- 
duction of a small concentration of vacancies,E3 certainly 
provided the average defect spacing is larger than extent 
of the most local ground-state degree of freedom in the 
pure system. Behaviour characteristic of the pure system 
may persist to much higher defect concentrations, since 
Dm/N > as long as more than three quarters of all 
sites are occupied. 
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